Classical noise and flux: the limits of multi-state atom lasers 
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By direct comparison between experiment and theory, we show how the classical noise on a multi- 
state atom laser beam increases with increasing flux. The trade off between classical noise and flux 
is an important consideration in precision interferometric measurement. We use periodic lOps radio- 
frequency pulses to couple atoms out of an _F = 2 *^Rb Bose-Einstein condensate. The resulting 
atom laser beam has suprising structure which is explained using three dimensional simulations of 
the five state Gross-Pitaevskii equations. 
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It is the high flux, spectral density, and associated first 
order coherence that has made the optical laser central 
to many technologies. In the field of precision measure- 
ment, atom lasers hold similar promise JJ. In a Sagnac 
interferometer, for example, the inherent sensitivity of a 
matter wave gyroscope exceeds that of a photon gyro- 
scope with the same particle flux and area by 11 orders 
of magnitude 2]. In any practical application of inter- 
ferometry to high precision measurement, whether it be 
with photons or atoms, there will be a trade off between 
the classical noise and the quantum noise or flux of the 
source. In this Letter, we investigate this trade off for an 
atom laser. 

In an unpumped laser, classical noise is the presence 
of unwanted excited dynamic modes. We find agreement 
between our experimental results and a full 3D Gross- 
Pitaevskii (GP) model, and show that at high flux, clas- 
sical noise increases with increasing flux. That we can 
achieve agreement between a 3D theory including all Zee- 
man states and experiment is significant. It is highly 
likely that the much sought-after pumped atom laser 
will operate under rather specific conditions of scatter- 
ing length, temperature and number density and ex- 
periments will need to be guided by accurate theoretical 
models that must be validated against experiments if we 
are to trust their detailed predictions. 

Mewes et al. demonstrated the first atom laser 
based on the application of pulsed radio- frequency (RF) 
fields to induce controlled spin flips from magnetically 
trapped to un-trapped states of a Bose-Einstein conden- 
sate. Later it was shown by Hagley et al. that a pulsed 
Raman out-coupling could be used to achieve a quasi- 
continuous multi-state atomic beam [^. Bloch et al. 
achieved continuous RF out-coulping for up to 100 ms, 
producing a single state atom laser beam, and showed 
that this beam could be coherently manipulated in di- 
rect analogy to the optical laser 0, Q • Both temporal 
and spatial coherence have been measured, and it has 
been demonstrated that RF outcoupling preserves the 
coherence of the condensate B 0, . The beam diver- 
gence has been measured [ll|, and there has been one 
real time measurement of the flux of an atom laser beam 



PJ. Ballagh et al. h^ introduced the Gross-Pitaevskii 
equation as an effective tool for investigating the atom 
laser within the semi-classical mean-field approximation 
and a number of groups foun d g ood agreement between 
GP models and experiment 0, Rabi cycling be- 

tween Zeeman components, a manifestation of the non- 
Mar kovian nature of the atom laser p^ll^ . was observed 
in the experiment of Mewes et al. and could be expected 
to significantly increase the amplitude, and possibly fre- 
quency, noise of the beam. There has been no investiga- 
tion of the relationship between classical noise and flux 
in an atom laser, and it is this aspect that we investigate 
both experimentally and theoretically in this Letter. 

Experimentally, a continuous atom laser based on res- 
onant output coupling put^ stringent limits on the sta- 
bility of cold atom traps (sj. Typical condensates have 
a resonant width of 10 kHz. Output coupling requires a 
stable magnetic bias, Bq, at the 0.1 mG level, one to two 
orders of magnitude better than typical magnetic traps. 
In comparison, a pulsed atom laser is relatively straight- 
forward to implement. A 10/xs pulse has a frequency 
width of 200 kHz, signiflcantly broader than both the 10 
kHz resonant width of the condensate, and the instability 
of our trap which fluctuates within the range ±15 kHz 
due mainly to thermal fluctuations of the coils. In the 
work reported in this Letter, we have opted to study a 
pulsed atom laser to ensure shot to shot reproducibility 
and allow detailed quantitative comparison to numerical 
models. We have chosen to study atom laser beams de- 
rived from an = 2, m^? = 2 condensate because of the 
richness and complexity offered by the system, although 
non-Markovian effects are present even in two-component 
atom lasers due to the nonlinear atomic dispersion rela- 
tions. 

In our experiment, we produce an F = 2,mF = '2 *^Rb 
condensate, consisting of approximately 50,000 atoms, 
via evaporation in a water-cooled QUIC magnetic trap 
13 with a radial trapping frequency Vr — 253 Hz, an 
axial trapping frequency i^z — 20 Hz, and a bias field 
Bq ~ 1 G. After evaporative cooling, the BEG was left 
to equilibriate for 100 ms. We then triggered an RF 
signal generator set in gated burst mode. The RF pulses 



2 



were amplified (35 dB) and radiated perpendicular to the 
magnetic bias field of the trap through a 22 mm radius 
single loop, approximately 18 mm from the BEC. To en- 
sure that we had calibrated all experimental parameters 
correctly, we made an initial series of measurements of 
the number of trapped and un-trapped atoms after the 
application of single RF pulses of varying amplitude; it 
was critical to establish agreement between experiment 
and theory in a simple mode of operation before pursuing 
studies of more complex dynamics. In Fig. 1, we show the 
results of these measurements in comparison with a one 
dimensional Gross-Pitaevskii theory of the atom laser de- 
rived from the full 3D model described by the following 
equations 

Z02 = {C + VT + Gy- 2A)02 + 2f)<?!>i 

i^i = (£ + iv^T + Gy - A)0i + 2VL(t)2 + V6f^0o 
i</>o ^{C + Gy)(j)o + V6n<j)i + V6n<j)-i 

= (£ - ^Vt + Gy + A)(/)_i + 2f}</>_2 + V617(/)o 
i0_2 = {C-VT + Gy + 2A)0_2 + 2r!0_i , 

(1) 

where is the GP function for the ith Zeeman state. 
Vt = \\x^ + ?/) + -iv2 + 

Here A — v^jv^ — 12.65 is the ratio of trapping frequen- 
cies, A and f2 are respectively the detuning of the RF 
field from resonance and the Rabi frequency, measured in 
units of LOz = 27rz^2 . U is the two-body interaction coeffi- 
cient and G =^ zamg / {fiu z) , where m is the atomic mass 
and g the acceleration due to gravity, zq = ^/hJnujJl 
is the usual harmonic oscillator length. The wave func- 
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FIG. 1: Outcoupled fraction as a function of RF amplitude. 
The solid square is the power at which the results of Fig. |21 
were generated. Theoretical curves are: Solid line, ttif = 2, 1 
trapped states, dot-dashed mp = and dashed, mp — —2, —1 
anti-trapped states. The experimental results are: triangles, 
niF ~ 2, 1 trapped states, circles mp = 0. Typical error 
bars are ±5% vertically and ±10% horizontally. Simulation 
parameter: A = —640. 



tions, time, spatial coordinates, and interaction strengths 
are measured in the units of zq, and uj~^. The excel- 
lent agreement between experiment and theory shown in 
Fig. 1 indicates that we have a well calibrated and repeat- 
able experiment. Up to the experimental uncertainty in 
the detuning A, there are no free parameters in the ID 
GP model. The theoretical results of the Rabi oscilla- 
tions presented in Fig. I are in good agreement with the 
approximate analytic theory presented by Graham and 
Walls for the limit of strong out-coupling . 

In the experimental data shown in Fig. 2, we present 
five pulse trains outcoupled from separate F=2, nip ~ 2 
condensates. In each case, the pulse train has been out- 
coupled in an 8 ms time frame. We wait 2 ms after the 
pulse window before turning the trap off to allow the final 
atomic pulse to completely separate from the condensate. 
After a further 2 ms, to allow expansion of the conden- 
sate, we image the condensate and pulses with a single 
lens onto a 12 bit CCD camera. For one, two, three and 
four RF pulses, we observe predictable out-coupling from 
the atom laser system. Figure 2(a) is indicative of this 
behavior, where four RF pulses (separation 2 ms) have 
been applied to the BEC, and we see four mp = atomic 
pulses in the positions expected from gravity. In 2(b) five 
RF pulses (separation 1.6 ms) have been applied and we 
observe five atomic wave-packets, again in the expected 
positions. However, we note that in the later 3 pulses 
there is a significant blurring with atoms appearing be- 
tween the expected positions of the pulses. This effect is 
not due to interference between the wavepackets. 

The transition from constant pulse amplitude shown 
in Fig. 2(a) to varying or noisy pulse amplitude with in- 
creasing repetition rate shows a clear trade off between 
classical noise and flux in the atom laser output. In 
Fig. 2(c) six RF pulses were apphed (separation 1.2 ms), 
however only five atomic pulses were observed, with the 
first atomic pulse being entirely absent. This observation 
is quite repeatable. The complete absence of a pulse is an 
extreme example of the trade off between classical noise 
and flux, and it is the dynamics behind this phenomenon 
that we wished to understand by comparison with a com- 
plete 3D GP model. At the higher pulse repetition rate 
used in Fig. 2(d) where the separation between pulses 
is 1 ms, the output is further distorted from the ideal. 
In Fig. 2(e), where the time between pulses has been 
reduced to 800/iS, the atom laser beam is longer than 
expected from pure gravitational acceleration. This can 
be explained by the influence of the the anti-trapped m p 
states on the mp = atoms that comprise the outcou- 
pled beam. It is quite clear from the data that increasing 
flux (and therefore decreasing shot noise) comes at the 
price of increasing classical noise. 

We have quantitatively modelled the experiment with 
a full 3D GP simulation including all five Zeeman states 
and with only the detuning A as a free parameter (Eq.^. 
This is a unique feature of the work presented here and 
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FIG. 2: (color) A series of pulsed atom lasers at diflFerent pulse 
rates. The applied radio-frequency (RF) pulses are varied 
from (a) 4 pulses, (b) 5 pulses, (c) 6 pulses, (d) 7 pulses, (e) 
10 pulses in an 8 ms window. The lower plots show a cross- 
section down the centre of the absorption data. The three 
dashed lines correspond in descending order to the centre of 
the condensate, the half-width of the 200 kHz RF resonance 
(100 kHz « 12/i m in our trap) and the position coinciding 
with the final RF out-coupling pulse (4 ms prior to imaging). 

allows us to understand all aspects of the GP physics, 
and hence the experiment. A ID model accurately de- 
scribes a single out-coupling pulse because it is essen- 
tially independent of the spatial structure. However we 
found that a full 3D simulation was needed to accurately 
track the spatio-temporal dynamics of a multi-pulse ex- 
periment. We simulated up to 3.2 ms, allowing three 
pulses for each case. Parallelised code was run on twelve 
processors of the APAC National Facility , requiring 
up to 800 hours of processor time per simulation. The 
numerical method was the psuedo-spectral method with 
Runge-Kutta split time step developed at the University 
of Otago li^l ■ Spatial grid sizes and time steps were mon- 
itored throughout the simulations to ensure the accuracy 
of the numerical solutions; e.g. the preservation of the 
normalisation. Spatial grids were grown, and time steps 
decreased, as required. At the end of simulations, spatial 
grids in the direction of gravity were 2048 points for the 
trapped (niF ~ 2,1) and anti-trapped (mp = —2,-1) 
Zeeman states, and 4096 for the atom laser output state 
{nip = 0). In the tight and loose transverse directions, 
128 and 32 points were used respectively. The corre- 
sponding spatial lengths were chosen so that both the 
momentum space and real space GP functions fit the 




FIG. 3: (color) Numerical simulations of the cases of Fig. 
2(c) and 2(b). Case 2(c) is to the left of vertical line. Each 
image shows the GP wavefunction density, in arbitrary units, 
integrated through the tight trap direction. Density is indi- 
cated by the colobar at the bottom of the figure, with black 
and white corresponding to zero and high density respectively. 
Each image is 120/im in both directions. From top to bottom 
the rows are: t — 1.2 ms, just before the 2nd RF pulse; just 
after the 2nd RF pulse; t — 2.4 ms, just before the 3rd RF 
pulse; just after the 3rd RF pulse, t = is the beginning 
of the first RF pulse. The columns from left to right are 
the niF = 2, 1, 0, —1, —2 states. The top-left image therefore 
shows the trapped position of the initial mp = 2 condensate. 
Case 2(b) is to the right of vertical line. Only the mp = 
state is shown. The rows are as before but the RF pulses occur 
a.t t — 1.6 ms and t — 3.2 ms. Each image is 120/im verti- 
cally and 140/^m horizontally. The rightmost image shows the 
mp — state density on a slice plane through the tight trap 
direction, just after the 3rd RF pulse. It is 40/xm horizon- 
tally and about 150/im vertically, allowing the first pulse to 
be seen. Simulation parameters: A — —633, Q = 457. See 
text for discussion. 

grid. This was about 40 fj,m and 140 fj,m. in the tight 
and loose trap directions, respectively. In the direction 
of gravity, it was 120 /Ltm for the mp ^ states, and 
twice that for the mp = state. Absorbing boundaries 
were used for the mp = 0,-1, —2 states. 

The simulations reveal that all five Zeeman states are 
involved in determining the final form of the atom laser 
output. The mp = —2 state has the least effect, as it 
is not strongly populated, and it quickly disperses in its 
anti-trapping potential. However the mp = — 1 anti- 
trapped state is highly populated and is directly involved 
in the loss of the initial pulse for the case of Fig. 2(c). 
The simulation in the left section of Fig. |31 shows how 
the first mp = atom laser pulse is destroyed by the 
second RF pulse: it transfers nearly all of the mp = 
component, produced by the first pulse, into the other 
four Zeeman states (second row, Fig. O). A new mp = 
pulse, somewhat lower than the first, originates from 
the mp = —1 state. However, it retains the momentum 
spread due to the anti-trapping potential, which causes 
it to disperse and be lost, so that it is not observed in 
the experiment. 
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The second atom laser pulse is in fact two distinct 
pulses; an upper one originating from the mp — 2 state, 
and a lower one from the mp — 1 state. This can be 
seen most clearly in the third row of Fig. O after they 
have become well separated. Since the two pulse com- 
ponents are not resolvable in the experiment, this is an 
example of the dynamics revealed by simulation. These 
components have different initial momenta. The mp = 1 
component, which originated from the mp — 2 state in 
the first RF pulse, was moving down towards its trap 
equilibrium when the second RF pulse arrived. 

The lower mp = Q pulse created by the second RF 
pulse escapes the fate of the first pulse because its down- 
ward momentum takes it lower than the first pulse, away 
from resonance. In fact, its position is close to that of the 
upper second mp = Q pulse during the third RF pulse in 
the case of Fig. 2(b), and it survives for similar reasons. 
This can be seen by comparing the bottom rows of the 
mp = Q columns of Fig. |3 Similarly the simulations ex- 
plain the relative intensity of the first and second mp = Q 
pulses in the experimental case of Fig. 2(b). 

We have shown, for the case of the F = 2 atom laser, 
that beyond a critical flux the classical noise on the out- 
put beam increases with increasing flux. The prospect of 
combining atom lasers with atom chips opens up enor- 
mous possibilities in precision measurement. Consider- 
ations of the trade off between classical noise and flux 
in atom lasers will be important in many applications in 
this field. We would expect many of the effects described 
here to be smaller for the F = 1 atom laser but not ab- 
sent. Rather than the two trapped states present in the 
F = 2 laser, only the mp = —1 state is trapped. Al- 
though atoms in the mp ^ 1 state are anti-trapped for 
the F = 1 laser, this state would be significantly pop- 
ulated for strong outcoupling and could be expected to 
contribute to classical noise on the mp — output beam 
just as the anti-trapped mp = —1 state does for the 
F ~ 2 laser studied in this Letter. The effects that we 
have described will be important not only for pulsed atom 
lasers, but also for unpumpcd continuous atom lasers. 
Just as in the pulsed case, at high flux, atoms will not 
only be coupled to the output beam, they will also be 
coupled to other trapped and untrapped Zeeman states 
and can be backcoupled from the output beam to the con- 
densate. The situation is complex and requires detailed 
investigation. The quantitative comparison between the- 



ory and experiment presented here is unique and points 
the way to the future development of atom laser sources 
for precision measurement. This is particularly true for 
the development of the pumped atom laser, one of the 
most important and sought-after devices in the field of 
atom optics. 
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